Ladder approximation to spin velocities in quantum wires 
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The spin sector of charge-spin separated single mode quantum wires is studied, accounting for 
realistic microscopic electron-electron interactions. We utilize the ladder approximation (LA) to 
the interaction vertex and exploit thermodynamic relations to obtain spin velocities. Down to not 
too small carrier densities our results compare well with existing quantum Monte-Carlo (QMC) 
data. Analyzing second order diagrams we identify logarithmically divergent contributions as cru- 
cial which the LA includes but which are missed, for example, by the self-consistent Hartree-Fock 
approximation. Contrary to other approximations the LA yields a non-trivial spin conductance. 
Its considerably smaller computational effort compared to numerically exact methods, such as the 
QMC method, enables us to study overall dependences on interaction parameters. We identify the 
short distance part of the interaction to govern spin sector properties. 

PACS numbers: 71.10.Pm,73.21.Hb,75.75.+a 
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At low energies the one-dimensional (ID) electron liq- 
uid exhibits charge-spin separationi*^, in contrast to 
the Fermi liquids of higher dimensionality^ Theoreti- 
cal considerations based on the Tomonaga-Luttinger liq- 
uid (TLL)^ have predicted this behavioni and recent 
experimenli^ has provided strong evidence supporting 
different velocities Vp and Va- Both velocities differ from 
the Fermi velocity vp: in the Hubbard model, for exam- 
ple, charge density waves (plasmons) propagate faster, 
Vp > vf, while spin density waves propagate slower, 
Va- < vf, than the Fermi velocity for repulsive on-site 
interactioni. 

Quantitative knowledge of is decisive to predict 
spin transport®'^ and magnetic propertiesiS. Within the 
random phase approximation (RPA), or when treating 
left and right going particles as distinguishableSiii as 
originally done by Luttingei^, Va = vf stays unrenor- 
malized with interactiongi^. On the other hand, even 
in first order perturbation theory the exchange or Fock 
term, proportional to the 2fcF-Fourier component of the 
electron-electron interaction V{q), influences the spin ve- 
locity. (In contemporary quantum wires the range of 
the interaction usually exceeds the Fermi wavelength, so 
V{q = 0) 7^ V{q — 2A;f)). Quantitatively, however, this 
latter approach is limited to V{2kF) <Si vf (correspond- 
ing to electron gas parameters rg <C 1). Furthermore, 
it spoils SU(2) spin rotation invariance of the micro- 
scopic electron model, see below. Quantum Monte Carlo 
(QMC) has demonstrated that v^/vf decreases with in- 
creasing inter-particle repulsion in quantum wires^'^, in 
qualitative resemblance to the Hubbard model^. Values 
of v„/vf ~ 0.5 have been estimated-?, for present day 
single channel quantum wiresi^ii^*i^. 

Numerically exact techniques, however, are computa- 
tionally extremely demanding, especially in spin sector. 
This yet has prohibited scans through larger parts of the 
parameter space that characterizes the microscopic in- 



teraction; in Ref. [13 only one interaction range and one 
channel width has been investigated. A future study of 
multi-channel quantum wires would require using con- 
siderably bigger systems. Here, a sufficiently accurate 
and tractable approximative scheme would be helpful. 
Among the techniques established for Fermi liquids the 
Hartree-Fock approximation (HF) has proven as useful 
to estimate boundary exponentaii and, when carried out 
self-consistently, to yield amazingly quantitative results 
in charge sector—; it captures for example very well 
the non-monotonous dependence of Kp on the electron 
densitjsiSi, beyond the RPA. On the other hand, the mean 
field approximation has turned out to fail badly in spin 
sectorplS.. Below we show that in a perturbative language 
this failure can be traced back to the wrong class of dia- 
grams summed by the self consistent HF. In the present 
work we sum up ladder diagrams and demonstrate that 
they comprise a 'complementary' class of diagrams that 
account much better for spin properties. 

The ladder approximation (LA) to the effective inter- 
action vertex originally has been established to study 
strongly correlated fermion systems with short range in- 
teractions, such as nuclear matter— but it performs re- 
markably well for the the also strongly correlated charge 
sector of a f DES2J-2a. 

In comparison to other Fermi- liquid methods-^, such as 
the Singwi, Tosi, Land, and Solander (STLS) approxima- 
tion schemeSiSiiS^, the LA is known to account well for 
short distance properties of the interaction when tested— 
with exactly solvable models2^. As we shall demonstrate, 
this short distance behavior is indeed most relevant to 
magnetic properties. In this paper we generalize the LA 
to allow for non-zero magnetizations and for spin cur- 
rents to gain homogeneous and static spin susceptibili- 
ties. Summing ladder diagrams yields estimates to 
that compare well with existing QMC-data at not too low 
electron densities while the numerical effort stays com- 



2 



parable to self consistent HF calculations. This enables 
us to scan different parameters of the microscopic inter- 
action. Of prime interest, e.g. to carbon nanotubes are 
the interaction range and the diameter of the quantum 
wires. 



I. MODEL 

To model the microscopic interaction we useiS 
2vp 



V{q) 



kpaB 



KoM)-Ko(gVrf2+4i?2) (1) 



in momentum space, where Kq denotes a modified Bessel 
function and ae is the Bohr radius. This form accounts 
for the experimentally important parameters, the inter- 
action range R, given by the distance to the nearest metal 
gate, and the diameter d of the wire, which determine 
V{q) at small and at large q, respectively. Details of the 
wire's cross section affect the interaction only at q > d~^, 
i.e. at carrier densities where occupation of higher sub- 
bands starts. We do not consider this case here. 

In Fermion representation the quantum wire (length 
L, s = ±1 denotes spin) is described by 



H 



- E 



''k — q^s 
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Ck,s are Fermi operators in the wave number basis. For 
non-linear single particle dispersion e(fc) = k"^ /2m, as in 
semiconducting quantum wires of effective carrier mass m 
and in the presence of 2fcF-scattering between antiparallel 
spins near opposite Fermi-points bosonization of model 
introduces terms of higher than quadratic order in 
the Bosonic density fluctuations. One could attempt, at 
least in principle, to eliminate higher order terms by the 
RG method in course of which they would renormalize all 
pre- factors of the quadratic terms, i.e. the TLL velocities 
VNf, and vji, in 



H. 



TLL 



q=iO 



L^ — (T,p 



4^ E i^N.N^ + vj,J^) 



(3) 

In practice such a RG approach to the quantitative val- 
ues of the TLL parameters has not been tested and seems 
not promising for strong interactions. On the other hand, 
there is little doubt that the one-dimensional electron 
liquid remains in the TLL phase even at strongest in- 
teractions. Fortunately, exact thermodynamic relations 
allow quantitative determination of the TLL velocities 
vnu and vj^, related with the total charge or spin N^, 
and the total currents J,y, from homogeneous and static 
susceptibilities since the latter are observable quantities. 

Symmetries reduce the number of independent TLL- 
parameters to be determined for given microscopic single 



particle dispersion e(fc) and interaction V{q). Left-right 
symmetry, e(fc) = e(— fc) leads to the TLL-rclations 



(4) 



for the sound velocities Vi, at which Bosonic charge or 
spin density fluctuations move, as described by the op- 
erators hi, q in Most quantum wires, furthermore, 
show Galilei invariance in charge sector which leads to 
vjp = v-p, independent of the interaction. Spin rota- 
tion SU(2) invariance, present in (0) without a Zeeman 
field or spin-orbit couplingSi, ensures the three velocities 
Va — VNa — vju to be cqual in spin sectoi^S. 

Here we concentrate on spin sector, where the thermo- 
dynamic relations read (see e.g. Ref. UXt) 



VNo 



and vjcr = 
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P = TrN„/2L and / = TrJg-/2L represent magnetization 
and spin current, respectively. 

For later comparison with the Hartree-Fock approxi- 
mation it is useful to express ground state energies and 
susceptibilities in terms of the self-energy 



i:,(fc,c^) = G°(fc,c^)-i -G,(fc,c^)- 



(6) 



Gs{k,uj) is the full 
tor with respect to (O 
{uj ~ k^/2m + sign[(fc - sl)^ - {kp + sPf]iO) 



electron propaga- 
and G°{k,uj) 

the 



free propagator. The ground state energy, required in 
(O, can be expressed aa^ 



(7) 



which reduces to 
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for self energies not depending on frequency lu (and there- 
fore are real as a consequence of the Lehmann representa- 
tion of Gs). This condition is fulfilled within the HF (c/. 
Eq. (|13|l below) and within the LA. We are not aiming 
to determine temporal or spatial correlation functions in 
the Fermionic representation, but rely here on the exact 
solution^ ^■2'' of the Boson model 

It can serve as a test to the quality of any approxi- 
mative scheme in spin sector whether or not it respects 
the conditioHpiS, v^a — vja dictated by SU(2) invariance. 
The HF, for example, leaves vj^ — f f independent of in- 
teractions since it ignores 2fcF-interactions between elec- 
trons of opposite spins and thus violates spin rotation 
symmetry. The STLS only allows to calculate magnetic 
susceptibilities v^a in spin sectoB^, so it cannot be tested 



3 



in its behavior regarding the SU(2) symmetry. The LA fi- 
nally does indeed rcnormalize ujo-, though only by about 
half the amount it renormalizes the value for v^a com- 
pared to vp. In so far we find that the LA does not fully 
obey the SU(2) symmetry but proves as superior to the 
other approximative schemes. 

II. LADDER APPROXIMATION 

To introduce the LA we follow the References and 
I21I but generalize the calculations for finite magnetization 
P and spin current /. The self-energy shift 

kp+s'il+P) 

Ss{k;P,I,kp) = J dk'x (9) 

-kF + s'{I-P) 

(^gss'{k,k\0)-5ss'gss'{k,k',k~k')^ 

due to interactions with the sea of other electrons is ex- 
pressed in terms of an effective (Brueckner) interaction 
matrix gss'{k, k' , q). Knowing gss'{k, k' , q) exactly would 
yield the exact ground state energy through and, by 
virtue of ((SJl, the exact values for the TLL- velocities. The 
Brueckner interaction matrix is closely related with the 
static structure factor which often is exploited to inves- 
tigate how the short range part of the interaction affects 
short distance correlations of one-dimensional electron 
liquids2L2a. 

For the ladder approximation a Bethe-Salpeter inte- 
gral equation 

00 

gss'{k,k',q) = ^(9) + ^ j <^P gss'{k,k' ,p) 

— 00 

xKss'{k,k',p)V{q~p) (10) 

has to be solved, describing multiple scattering between 
electrons that otherwise propagate freely outside the 
Fermi sphere according to 

K,,,{k,k',p)= (11) 

^ e{\k~si^p\-kFs)e{\k' ~s'i + p\-kFs') 
k^ + k'^ - {k - p)^ - {k' + p)^ 

In (|12ll we have defined spin dependent Fermi-momenta 
kps = kp + sP at finite magnetizations. We solve equa- 
tion H10|l for each pair of momenta (fc, k') and spins 
(s, s') the Householder method after mapping the infi- 
nite p range to [0, tt] by p = kp cot ((> and discretizing 
= Tr{i - 1)/{N - 1). A value of TV = 150 turned out 
as an optimal compromise between CPU time (increas- 
ing with N^) and accuracy (the error in g ^ N^"^ due 
to the cusps of g{k,k',q) at g = 2fcF)- Only at strong 
interactions, when spin velocities drop below half of the 
Fermi velocity, we used N = 250. Finally, ground state 
energies are integrated via Q and © using the trapezoid 
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FIG. 1: Self-consistent g-fi(k,k',0) for R/d = 14.14, kpd = 
0.15, P = 0.3,/ = 0.2. Due to correlation eflfects the 
Brueckner interaction vertex shows pronounced valleys along 
k + k' = ±2P. The bare value V{0) = 22.29t;F. 




FIG. 2: Ground state energy £0 = Eo/(Lkp/m) versus fcpd 
for R/d = 14.14. The electron density 2k-e/'n can be ex- 
pressed through the electron gas parameter = 7r/(8fcFaB) 
at d = 'LA' are our results, 'pert' denotes first order 

perturbation theory when g(k,k',q) = V(q)- 'HF' refers to 
self-consistent Hartree-Fock results of Ref. [13 and the QMC 
data are taken from Ref. IT^ 

rule on a grid of 129 points over the interval [—2kF; ^kp] 
which sufficed to accurately resolve the effects of small 
magnetizations . 

The solution gss'{k,k' ,q) exhibits pronounced depen- 
dence on k and A:', cf. Figure ^ arising primarily from 
the short distance correlations at k = — fc'. Ignoring 
this dependence, as it is commonly done in 3D {cf. Ref. 
29 and references therein), would clearly not be justi- 
fied. This correlation has similarities to the striking 
anti-ferromagnetic 2fcF modulations found in the self- 
consistent HF-density in ID^^ which also cannot be ig- 
nored for reliable results in the charge sector. 

III. RESULTS 

Figure |21 shows ground state energies (in units of the 
respective kinetic energies) versus carrier density using 
different approximations. QMC datai^ can be regarded 
as exact within symbol size. Remarkably, at densities 
kpd ^ 0.2 , where self-consistence is seen in Figure |21 to 
still improve the HF-estimate to Eq compared to 'pert' 
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FIG. 3: Charge sector Luttinger parameter Kp versus fepd for 
R/d = 14.14. Labels as in Figure |21 
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FIG. 4: Spin velocities in units of up versus fcpd for R/d — 
14.14. The HF data are taken from Ref. [3 and the QMC 
data from Ref. 0. 

(and then provides the optimum Fermi-liquid wave func- 
tion), the LA yields slightly lower ground state energies. 
As seen in Figure O the LA also yields better values of 
Kp in this regime. Only at smaller electron densities HF- 
theory engenders lower ground state energies and a better 
estimate to the /Cp-parameter {kpd < 0.1). This energy 
gain is accompanied by spontaneous symmetry break- 
ing and pronounced (though unreal) static 4A:F-periodic 
Wigner crystal-like modulations^^ of the charge density 
in the HF ground state. However, this success of the HF 
in charge sector does clearly not carry over to the spin 
sector— as seen in Figure ^ 

Remarkably, in the LA the spin velocity follows the 
QMC-data down to rather low electron densities, supe- 
rior to other electron gas theory approaches. Only at 
small densities vanishes, pretending the transition 
into a ferromagnetic ground state which is not expected 
to occur in ID for finite range interactions. The LA 
does not reproduce the behavior Va- oc derived from 
the Hubbard model at small fillings'^" to which quantum 
wires should cross over— when the inter-particle spacing 
exceeds the interaction range, k-pR <^ 7r/2. 

Figure El demonstrates that at interaction ranges 
kpR 3> 7r/2 the value of R does not affect the spin 
velocity, unlike the ground state energy or the charge 
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FIG. 5: Spin velocity versus interaction range R for fcpd = 
0.3. At the arrow R = 7r/4fcF = 2rsaB. When R <ti d the 
interaction Q vanishes trivially; we always assume R ^ d. 
The inset shows spin velocity versus wire width for kpd = 0.3 
at R/d = U.U. 

sector exponent which both depend logarithmically on 
R/di2i&. This is consistent with the perturbative re- 
sult according to which V{k — 2kp) but not V{k — 0) 
governs the magnetic susceptibility where the former 
V{2kF) ^ 2vpKo{2kpd)/k-paB becomes independent of 
kpR 3> 1. This is of particular relevance to carbon 
nanotubea^ where kpR can be of the order of 10'^. The 
inset of Figure|Slcomplements Figure^ showing how the 
spin velocity varies with wire width d at fixed kpd — 0.3. 

IV. DISCUSSION 

To analyze further why the LA captures magnetic 
properties so well we compare it diagrammatically with 
the self-consistent HF. For this purpose we use the self- 
energy (jHJ were HF and LA can be compared directly. We 
expand the self-energy @ to second order in the inter- 
action. The first order contribution, obtained by putting 
g{k,k',q) = V{q) in ®, yields 

fcp+sP 

i:.,pc..t(fc + .s/) = ^F(o)- j dk'^^^^^. 

-kF-sP 

(12) 

The exchange term on the right hand side of H12|l is effec- 
tive only for parallel spins and therefore independent of 
the spin current / so that vj^ remains un-renormalized 
to this order. 

Self-consistent HF is described by 

^.,HF(fc) = ^Y.J d^'^^ ^ 

iV{0)-Sss'V{k-k'))Gs'{k',Lu) (13) 
Gs{k,uj) = G°(fc,c^)-HG°(fc,c^)i:,3F(fc)G,(fc,c^)(14) 

To first order, Us,piF{k) agrees with ()12(l . We see that 
(|13|l will not depend on /, even when self-consistence 
G" ^ G is reached. 
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FIG. 6: Irreducible electron self-energy S{k) in second order perturbation theory. The LA accounts for the contributions (a) 
and (fe) while the HF includes the terms (c) — (e). At fc = fcps the term (a) behaves like PlnP yielding a divergent contribution 

to Va-. 



Second order contributions are obtained after iterating respectively. For the LA this yields 
® and (Unjl for the LA or ^ and ^ for the HF, 



(a) + / Ak'dujG°ik',uj) J dpdcj' V{pfG°,{k ~ p, /2m ^ ij')G°,ik' + p, k'^ /2m + uj') 



(15) 



1 



(2^)4 
while for the HF 

'^s.Hf(^) ~ '^s,pcrt(fc) 



dk'duj G°(k', oj) / dpduj' V{p)V{k - fc' - p)G°(fc - p, k^ /2m - u')Gl{k' + p, k'^ /2m + J) 



(16) 



(d) 

(e) 
(/) 



j^,V{Of j dfc'd^G" (fc^',c.) j dpduj' Gl„{p,u:'f 
(^nO)E / dfc'dc.GO(fc',c.) j dpdLo'Vik' -p)Gl.{p,oj'f 



1 



{2^Y 



dk'dLoGl{k',uj) / dpdcj'V^(fc-p)V^(p-fc')G°(p,t^')^ 



(^nO)X^ / dfc'dc.G°(fc',c.) y dpdc.'y(fc-p)G°(p,c.')' 
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These are all irreducible second order contributions. 
None of the diagrams, depicted in Figure El occurs in 
both of the approximations, insofar the LA and the HF 
can be regarded as complementary. It is knownii that 
any finite (but the first) order self-energy contribution 
diverges logarithmically for fc ^ fcp ; to second order this 
occurs only due to the 'LA terms' (a) and (&) while the 
HF contributions remain finite. 

Furthermore, only term (a) breaks Galilei invariance 
in spin sector as required to satisfy the SU(2) symme- 
try condition according to which the spin conductance 
should be equally renormalized by interactions as the 
magnetic susceptibility, vja = Vno-- Term (a) in H15(l 
describes scattering between antiparallel s = — s' spins 

(2) 

at opposite Fermi points which affects djUs (fc; P, I, fcp) 
and the spin conductance vj^- 



In order to compare the importance of the different 
terms (a) to (/) regarding dpSs (fc;P, /, fcr) and thus 
the magnetic susceptibility VNa we consider for the mo- 
ment a contact interaction, V{q) — V. Then all second 
order contributions can be calculated analytically. The 
sum of all spin parallel parts s = s' cancel (spin parallel 
Fermions do not interact at contact) and only the s = —s' 
part of (a) and the term (c) (with s' — s or s' — s") re- 
main non- vanishing. The latter is not divergent while 
term (a) can be expressed by sums of dilogarithms2^ 
which can further be analyzed for I,P <^ fcp. The re- 
sult is 

(a) oc mV^[7r^/3 - i\P + I\)/kFH{\P + iD/kp) - 

2{P + I)/kF + 0{P^,I^)] (17) 

at fc = fcp + sP + si and / ^ — / in (|17|) near the other 
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Fermi point k = —kp — sP + si. The infinite slope seen 
in equation I^TJ at / = P = (by virtue of © and ©) 
results in logarithmically diverging second order contri- 
butions to vjcr — v-p and to VNa—vp as a function of P or / 
which dominate over the non-diverging HF-term (c). The 
LA comprises in summing the leading logarithmically di- 
verging contributions to the spin velocity which is one 
reason for its success regarding spin sector properties. A 
finite interaction range does not remove logarithmic di- 
vergencies to qualitatively alter this observation. 

V. OUTLOOK 

In conclusion, we have generalized the ladder approxi- 
mation (LA) to investigate the spin sector of single chan- 
nel quantum wires in the presence of a realistic micro- 
scopic interaction. While the numerical effort is consid- 
erably smaller, we obtain values for the spin velocities 
that compare well with existing quantum Monte-Carlo 
data at not too small particle densities. The LA accounts 
for interaction diagrams which renormalize the spin con- 
ductance. Furthermore, the LA diagrams include the 
leading logarithmically diverging contributions to spin 
conductance and susceptibility and thus are of dominant 
importance for magnetic properties. The self-consistent 
Hartree-Fock approximation (HF) misses these diagrams 



and therefore leaves the spin conductance unaffected by 
interactions leading to an erroneous result for the spin 
velocity. The short distance part of the interaction is 
identified to govern spin sector properties. Metal gates 
fabricated close to a quantum wire will screen the long 
range part of the interaction but leave spin properties 
almost unaffected. 

The LA should prove useful to study spin properties of 
multi-channel quantum wires where the numerically ex- 
act methods would be even more demanding than already 
in the single channel case. Compared to the HF the LA 
does not suffer from incommensurate densities fep-f 7^ k-pi 
at finite magnetizations. This advantage over the HF 
might even carry over to the charge sector of coupled 
quantum channels since it is likely that incommensurate 
particle densities in different channels will cause similar 
instabilities of the HF-procedure as for magnetic proper- 
ties in single channelsiS. Thus, the LA might prove as 
the method of choice also for the charge sector of multi- 
channel quantum wires. 
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